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ABSTRACT: Knowledge of the molecular frontier levels’ align¬ 
ment in the ground state can be used to predict the photocatalytic 
activity of an interface. The position of the adsorbate’s highest 
occupied molecular orbital (HOMO) levels relative to the sub¬ 
strate’s valence band maximum (VBM) in the interface describes 
the favorability of photogenerated hole transfer from the VBM 
to the adsorbed molecule. This is a key quantity for assessing 
and comparing H2O photooxidation activities on two prototypical 
photocatalytic UO2 surfaces: anatase (A)-TiO2(101) and mtile 
(R)-TiO2(110). Using the projected density of states (DOS) from 
state-of-the-art quasiparticle (QP) Go Wo calculations, we assess 
the relative photocatalytic activity of intact and dissociated H2O 
on coordinately unsaturated (Ti CM5 ) sites of idealized stoichiomet¬ 
ric A-TiO2(101)/R-TiO2(l 10) and bridging O vacancies (0^ c ) of defective A-TiO2-jc(101)/R-TiO2-jc(l 10) surfaces (x = Vi, Vs) for vari¬ 
ous coverages. Such a many-body treatment is necessary to correctly describe the anisotropic screening of electron-electron interactions 
at a photocatalytic interface, and hence obtain accurate interfacial level alignments. The more favorable ground state HOMO level align¬ 
ment for A-TiO2(101) may explain why the anatase polymorph shows higher photocatalytic activities than the rutile polymorph. Our 
results indicate that (1) hole trapping is more favored on A-TiO2(101) than R-Ti02(l 10) and (2) HO@Ti CM s is more photocatalytically 
active than intact H20@Ti CM5 . 
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1. INTRODUCTION 

UO2 is widelv used in photocatalysis, photoelectrocatalysis, and 
photo voltaic s.^ 6 In particular, the H20-Ti02 interface has been 
intensively studied both experimentally 7 and theoretically. ^This is 
due to both the ubiquitous nature of the aqueous environment, and 
the technological importance of water splitting. 9 10 Because large 
single-crystal samplesy)f the anatase polymorph are less stable than 
the rutile polymorph, 11 13 most surface science studies have fo¬ 
cused on the rutile (110) surface of Ti02- However, in the nanopar¬ 
ticle form the anatase polymorph is more stable, 14 and moreover it 
has a higher photocatalytic activity.^ 

A proper assessment of an interface’s photocatalytic activity re¬ 
quires an accurate description of its frontier levels’ alignment. This 
is because interfacial electron transfer is controlled by the align¬ 
ment of the highest occupied and lowest unoccupied molecular 
orbitals (HOMO/LUMO) relative to the valence band maximum 
(VBM) and conduction band minimum (CBM). 45 In particular, 


H2O photooxidation, i.e., the oxygen evolution reaction (OER), is 
initiated by the transfer of the photogenerated hole from the sub- 
strate’s VBM to the HOMoP® 

Based on the HOMO’S position relative to the VBM prior to ir¬ 
radiation, i.e., vertical alignment, one may establish trends in pho¬ 
tocatalytic activity among a group of systems. 19 20 Even in cases 
where the HOMO initially lies below the VBM, 21 after light ab¬ 
sorption and nuclear relaxation, these levels may reorder, with the 
hole localized on the molecule. 22 Essentially, the closer to the 
VBM and more localized on the molecule the HOMO is initially, 
the greater the molecule’s propensity for trapping the hole. For 
these reasons, the alignment of the H2O occupied levels prior to 
irradiation is most relevant for understanding the OER. 

Recently, we applied many-body quasiparticle (QP) GW tech- 
niques 1 — L to determine the H2O occupied levels’ alignment on 
rutile (R)-TiO2(110). 1 —'We found that an accurate description of 
the interfacial anisotropic screening via QP GW is essential to ac- 
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curately describe the interfacial level alignment.^ 1 25 27 Specifi¬ 
cally, the occupied QP density of states (DOS) projected onto the 
molecule is an effective means for interpreting difference spectra, 
i.e., the difference between spectra with a chemisorbed molecular 
overlayer and a clean T 1 O 2 surface, from ultraviolet photoemission 
spectroscopy (UPS). 1 — 1 Such theoretical approaches are necessary 
to disentangle highly hybridized adsorbate levels from those of the 
substrate, such as those of the H 20 -Ti 02 interface. 1 — 1 

Here, we investigate the H 2 O occupied levels’ alignment on the 
anatase (A)-TiO2(101) surface, as it is the most common surface in 
nano structured Ti02- 22829 In the absence of UPS measurements 
for H 2 O on A-TiO2(101), we compare the results to the Go Wo 
PDOS of H 2 O on R-TiO2(110),® which is consistent with UPS 
difference spectra. 31 33 

In particular, we pe rform GoWo 232434 and partially self- 
consistent 35 (sc)QPGWl^- 26 calculations based on Kohn-Sham 
(KS) levels from density functional theory (DFT) using a local den- 
sity approximation (LDA ), 1 — 1 generalized gradient approximation 
(PBE/“ or a range-separated hybrid (HSE)® 39 exchange corre¬ 
lation (xc)-functionals. From these calculations we obtain the to¬ 
tal and projected QP DOS for a variety of coverages [Va to U /2 
monolayer (ML)] of intact and dissociated H 2 O adsorbed on co- 
ordinately unsaturated Ti sites (H20@Ti CM5 ) of stoichiometric A- 
TiO2(101) and bridging O vacancies (H20@0^ c ) of defective A- 
TiO2_i/ 4 (101) and A-TiO2_i/ 8 (101) surfaces with AML and 'AML 
O v b a r c . The Ticus and 0^ r sites of A-TiO2(101) and R-Ti02(l 10) are 
shown schematically in Figure [l] 


(a) A-TiO 2 (101) (b) R-TiO 2 (110) 



Figure 1 . Schematics of clean stoichiometric (a) A-TiC>2(101) and (b) R- 
TiC>2(110) surfaces. Ti and O atoms are depicted in silver and red, respec¬ 
tively. Coordinately unsaturated Ti sites (Ti CM5 ) and bridging O atoms (Otr) 
are labelled. 


On the one hand, by considering the absolute interfacial level 
alignment, i.e., relative to the vacuum level E vac , one obtains the 
interface’s ionization potential IP = -£vbm + ^vac- This is the quan¬ 
tity that can be compared with red-ox potentials.^ 041 Moreover, 
the absolute level alignment allows a direct comparison between 
alignments across different substrates,® such as A-TiO2(101) and 
R-TiO2(110). Finally, from the absolute level alignment, one can 
determine whether changes in the H 2 O level alignment across dif¬ 
ferent substrates are attributable to shifts of the substrate or molec¬ 
ular levels. 

On the other hand, by considering the level alignment relative 
to the YBM of the substrate £y bm, one may directly compare the 
favorability of photogenerated hole transfer from the substrate’s 
VBM to the molecule’s HOMO. In combination with the IP this 
allows a robust comparison of photocatalytic activity across sub¬ 


strates. Moreover, by referencing the spectra to the VBM, one may 
directly compare the shape and dispersion of the valence band edge 
at the VBM. Finally, the VBM is the most reliable KS energy ref¬ 
erence, from a theoretical perspective.®^ 

For these reasons, we shall make use of both VBM and E vac en¬ 
ergy references as appropriate. In particular, we provide the abso¬ 
lute level alignment when comparing to HSE DFT and HSE Go Wo 
DOS. This is because the VBM and CBM from HSE DFT for 
clean 26 and 1ML H20@Ti CM J 21 * 40 " 41 on R-Ti02(l 10) are consis¬ 
tent with measurements for clean quasi-stoichiometric 26 42 45 and 
liquid H 2 O covered R-TiO2(110), 41 respectively. 

We begin by providing a detailed description of the techniques, 
parameters, and terminology employed throughout in Section [2] 
In Section [3T] we compare our results to calculated and measured 
electronic band gaps E g , optical band gaps h(o g , and macroscopic 
dielectric constants £00 of bulk A-Ti 02 and R-Ti 02 available in the 
literature. To further test the robustness of our approach, and the pa¬ 
rameters employed, we compare the dielectric function obtained by 
solving the Bethe-Salpeter equation (BSE) based on Go Wo eigen¬ 
values with that obtained from reflection spectra measurements for 
bulk A-Ti02- We also compare the IP for clean idealized stoichio¬ 
metric A-TiO2(101) and R-TiO2(110). In Section |T2| we provide 
a detailed comparison of the H 2 O PDOS for intact and A dissoci¬ 
ated H20@Ti CM 5 on A-TiO2(101) and R-Ti02(l 10) relative to E vac 
and £ V bm with PBE DFT ’ HSE DFT > PBE scQPGWl, PBE G 0 W 0 
and HSE GqWq. In Section [33] we provide a similar detailed com¬ 
parison for dissociatively adsorbed H20@0^ c on A-TiO2-i/ 4 (101) 
and R-TiO2_i/ 4 (110) relative to E vac . Finally, in Section [T4| we ex¬ 
tend the comparison of A-TiO2(101) and R-TiO2(110) to their ad¬ 
sorption energies and level alignments relative to £vbm with H 2 O 
coverage (Va to U/ 2 ML), H 2 O dissociation (intact to fully dissoci¬ 
ated), and surface composition (0^ c coverage). This is followed 
by concluding remarks. 

2. METHODOLOGY 

2.1. Theoretical Methods. DFT based on standard xc- 
functionals, e.g., LDA and PBE, tends to significantly underes¬ 
timate the electronic band gaps of semiconducting metal oxides, 
such as TiC> 2 . 46 This is in part due to their underestimation of the 
screening of the electron-electron interaction. 

DFT based on hybrid xc-functionals, e.g., HSE, partially reme¬ 
dies this by replacing a fraction of the exchange term with Hartree- 
Fock exact-exchange. The fraction of Hartree-Fock exact-exchange 
included, a , acts as an effective constant screening of the Hartree- 
Fock electron-electron interaction term, i.e., an inverse dielectric 
constant for the system a ~ s ~^>. 47 In particular, we use the HSE06 
variant, with a range separation parameter of yu = 0.2 A -1 , of the 
HSE hybrid xc-functional, which includes 25% exact-exchange 

[To! 

(a = 0.25). 1 — 1 For systems where the screening is rather homoge¬ 
neous, and £00 ~ 4, this leads to a better description of the electronic 
band gap, 47 e.g., of bulk Ti 02 - 46 

However, since HSE applies the same screening to all the lev¬ 
els regardless of their nature, it fails to describe the anisotropic 

screening felt bv molecular levels at an interface. As a result, lo- 

I 21 I 25 I 26 I 

calized occupied molecular levels are underbound by HSE “— 1 — 1 
This can lead to significant errors in HSE’s description of molec¬ 
ular/semiconductor interfacial level alignment. 21 25 26 Instead, QP 
techniques, e.g., Go Wo and scQPGWl, where the spacial depen¬ 
dence of the screening is explicitly included, provide a better de- 
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scription of the interfacial level alignment. 1 —"— 1 

In the Go Wo approach, the contribution to the KS eigenvalues 
from the xc-potential V xc is replaced by the self energy E = /GW, 
where G is the Green’s function and W is the spatially depen- 
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dent screening 23 obtained from the KS wavefunctions.^ 4 The di¬ 
electric function is obtained from linear response time-dependent 
(TD) DFT within the random phase approximation (RPA), includ¬ 
ing local field effects.^ From GqWq one obtains first-order QP cor¬ 
rections to the KS eigenvalues, but retains the KS wavefunctions. 
Generally, these QP corrections to the occupied levels are linearly 
dependent on the fraction of the KS wavefunction’s density in the 
molecular layer. 1 1 ’ This means the screening of these levels W is 
quite anisotropic. For this reason, QP GW methods are necessary 
to accurately describe the interfacial level alignment. 

Moreover, to include dependencies on the QP wavefunctions, 
and possibly obtain an improved absolute level alignment for 

the interface, one can also emnlov self-consistent OP techniques, 

I21I25I26I 

such as scOPGWl. 1 1 1 J Here, we have employed a single-shot 
scQPGWl approach,^- 35 where 100% of the DFT xc-potential is 
replaced by the QP self energy in a single self-consistent scQP GW 
cycle. We employ this procedure to make practical scQPGW cal¬ 
culations for large interfaces, such as H2O-A-TiO2(101). In so do¬ 
ing, one obtains eigenvalues comparable to those from Go Wo, along 
with the QP wavefunctions. This differs from the scQPGWl ap- 
proach as previously applied to the H 2 O-R-TiO 2 ( 110 ) interface, 1 — 1 
where 25%, 25%, and 50%, of the QP self energies were “mixed” 
with the DFT xc-potential over three scQPGW cycles,^ respec¬ 
tively. 

2.2. Computational Details. Our G 0 W 0 calculations 23 24 34 
have been performed using vasp within the projector augmented 
wave (PAW) scheme. 48 The Go Wo calculations are based on KS 
wavefunctions and eigenenergies obtained from DFT using either 
LDA^PBE^or HSE^xc-functionalsP 

The geometries have been fully relaxed using LDA 3 ^or PBE, 37 
with all forces < 0.02 eV/A. The HSE 38 calculations have been 
performed for the relaxed geometries obtained with PBE. We have 
employed a plane-wave energy cutoff of 445 eV, an electronic tem¬ 
perature of kftT » 0.1 eV with all energies extrapolated to T -► 0 
K, and a PAW pseudopotential for Ti which includes the 3 s 2 and 

3/ semi-core levels. All calculations have been performed spin 

I2TI 

unpolarized. 1 — 1 

For the clean stoichiometric A-TiO2(101) surface we have used 
a five-layer slab, an orthorhombic lxl unit cell of 10.23 x 3.78 x 40 
A , a T-centered 4x8x1 k-point mesh, and approximately 9 l A un¬ 
occupied bands per atom. For the clean defective A-TiO2_i/ 4 (101) 
surface we have used a monoclinic 1x2 unit cell of 10.23 x 7.56 x 40 
A and a T-centered 4x4x1 k-point mesh. For the clean de¬ 
fective A-TiO2_i/ 8 (101) surface we have used a 1x4 unit cell of 
10.23 x 15.13 x 40 A 3 and a T-centered 4x2x1 k-point mesh. 
For the H 2 O covered surfaces, we have employed a five-layer slab 
with adsorbates on both sides, an orthorhombic lxl unit cell of 
10.23 x 3.78 x 47 A 3 , a T centered 4x8x1 k-point mesh, and ap¬ 
proximately 9Ye unoccupied bands per atom, i.e., including all lev¬ 
els up to 30 eV above the VBM, an energy cutoff of 80 eV for 
the number of G-vectors, and a sampling of 80 frequency points 
for the RPA dielectric function. The Go Wo parameters are con¬ 
sistent with those previously used for d escrib ing bulk R-Ti02, R- 
TiO2(110) clean surface and interfaces. 25 26 Although our Go Wo 
calculations do not include electron-phonon 50 and lattice polariza¬ 
tion 51 contributions, these parameters have been shown to provide 
accurate descriptions of bulk optical absorption spectra, and both 
clean surface and interfacial level alignment. 25 26 

It has previously been shown 46 3 that the experimental op¬ 
tical spectra for bulk A-Ti02 may be obtained via BSE 54 based 
on Go Wo eigenvalues. In our BSE calculations, we include the 
electrostatic electron-hole interaction using the effective nonlocal 
frequency independent exchange correlation f xc (r,r ',oj = 0) kernel 
suggested in ref. [55] For bulk A-Ti02, we have used a tetrago¬ 


nal conventional 12 atom supercell with experimental lattice pa¬ 
rameters a = b = 3.78 A and c = 9.5 A, 56 and a dense T-centered 
10x10x4 k-point mesh. For bulk R-Ti02, we have used a tetrag¬ 
onal 6 atom primitive cell with experimental lattice parameters 
a = b = 4.5941 A and c = 2.958 A, 56 a T-centered 6 x 6 x 10 k-point 
mesh with PBE and HSE and a denser T-centered 8 x 8 x 12 k-point 
mesh with LDA. For both A-TiC >2 and R-Ti02, we have included 
ftunocc = 12 unoccupied bands per atom. For the BSE calculations 
of bulk A-Ti02, we have used 480 sampling points for the RPA 
dielectric function, and included all the transitions between the 16 
highest energy occupied bands and the 12 lowest energy unoccu¬ 
pied bands] 

2.3. Terminology. To compare the relative stabilities of the 
H 2 O covered anatase and rutile polymorphs, we have performed 
single-point rpbeP based DFT calculations using the PBE re¬ 
laxed structure for the H 2 O adsorption energies E a d s on stoichio¬ 
metric A-TiC>2(101) and defective A-TiC^-jcQOl) surfaces. The 
RPBE xc-functional was especially developed for the prediction of 
adsorption properties on metal surfaces. 57 ^Furthermore, RPBE has 
been shown to provide accurate formation energies for metal diox¬ 
ides 58 and perovskites. 59 

The H 2 O adsorption energy on the T\ CU s site of a stoichiometric 
A-TiO2(101) surface is given by 

Eads « £["H 2 Q + A-TiQ 2 (101)] -E[A-TiQ 2 (101)] 
n 

( 1 ) 

where n is the number of adsorbed H 2 O functional units in 
the supercell, and E[wH 2 0 +A-TiO 2 (101)], £[A-TiO 2 (101)], and 
£[H 20 ] are the total energies of the covered and clean stoichio¬ 
metric surfaces and gas phase water molecule, respectively. Sim¬ 
ilarly, the H 2 O adsorption energy on the O v b a r c site of a defective 
A-TiO2-*(101) surface is given by 


Eads * £[rcH 2 Q + A-TiQ 2 - x (101)] -£[A-TiQ 2 - x (101)] ^ Q] 

n 

( 2 ) 

where Zs[nH20 +A-TiO2-*(101)] and E[A-TiO2-jc(101)] are the 
total energies of the covered and clean defective surfaces, respec¬ 
tively. 

To provide a quantitative comparison between the DOS for the 
H 20 -A-Ti 02 and H 20 -R-Ti 02 interfaces, we employ the inter¬ 
faces’ IPs. These are obtained from the difference in energy be¬ 
tween the vacuum level Evac and the VBM £vbm> /P = -£vbm + 
Evac , where Evac is the maximum surface averaged electrostatic 
potential in the vacuum region between slabs. 

Similarly, to provide a quantitative comparison between the 
PDOS for the H20-A-Ti02 and H20-R-Ti02 interfaces, we em¬ 
ploy both the highest H 2 O PDOS peak s P ^ S an d the average 
energy of the highest energy electron, or HOMO, of the PDOS, 
£ homo* To obtain s p ^ s from the PDOS, we fit three Gaussians to 
the first few peaks below the VBM. In this way we may disentangle 
the highest energy peak when it forms a shoulder within the upper 
edge of the PDOS. 

However, to assess trends in the comparative photocatalytic ac¬ 
tivity of the H 20 -A-Ti 02 and H 20 -R-Ti 02 interfaces, one should 
consider not only a peak’s energy, but also differences in its’ in¬ 
tensity, i.e., localization on H 2 O. Both quantities are incorporated 
within the single descriptor ^homo- We define £ homo as th e first 
moment of the PDOS, p PD0S (s ) over the interval encompassing 
the highest energy electron. More precisely, 


PDOS 

£ homo 


r 

JEi 


6TVBM+A 


PDOS/ \ / 
sp ( s)ds , 


(3) 


where eyBM ^e VBM energy, A ~ 1 eV ensures the tail of the 
VBM is included within the integral, and E\ is the lower bound of 
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the energy range encompassing the highest energy electron of the 
PDOS, i.e., 

rs w bm+A pD0S 

/ p (s)ds= 1. (4) 

JE i 

3. RESULTS AND DISCUSSION 

3.1. Bulk and (101) Surface of Anatase Ti0 2 . To test the re¬ 
liability of the parameters we have employed to calculate the Go Wo 
levels of A-Ti 02 , we first consider the optical response of bulk 
anatase. Previous DFT band structure calculations 46 52 53 found A- 
Ti02 has an indirect electronic band gap between the VBM along 
the 2 path at 0.88r - M®, i.e., 2, and the CBM at C Our PBE 
Go Wo calculation yields an indirect electronic band gap for A-Ti 02 
of 3.86 eV, from a VBM at 0.8T -> M. This is comparable with the 
Go Wo indirect band gaps reported in the literature, as shown in Ta¬ 
ble [l ] 

Table 1. Direct and Indirect Band Gaps E g and Optical Gaps ha> g in 
eV of A-Ti02 and R-Ti02. 


method xc-functional 

A-Ti0 2 

R-Ti0 2 


electronic band 

r^r E^r r^r 

gap 

' T^R 

DFT HSE 

3.72 3.63 

3.40 

3.40 

3.60“ 

3.39“ 

3.39“ 

LDA 

G 0 Wq 

3.93 3.86 

3.33 

3.26 

4.14* 3.56* 

3.73“ 

3.38* 

3.46“ 

3.34* 

PBE 

4.29 c 3.83“ 

3.59“ 


PBE+A 

3.57 d 

3.30 d 

3.23 d 

PES/IPES 


3. 

.3 ±0.5“ 

(axM-ha )) 2 

3.53 / 


3.37 / 


optical gap 


- PBE 

3.73 


3.15 

3.57“ 


3.28“ 

Transmission 

3.42* 



Absorption 



3.03* 

Reflectance 

3.21'' 


3.00‘ 


a Rcf. [46] ^Ref. [53] c Ref. [52] ^Ref. |60] c Photoemission and Bremsstrahlung isochromat spec¬ 
troscopy from ref. 61 ^Estimate assuming a nearly direct band gap based on Kubelka-Munk 
adsorption coefficients qtkm from reflectance measurements of phase-pure nanoparticles in ref. 
[62] 8 Ref. [63] ^Refs. |64| and [65| 'For pure-phase nanoparticles from ref. |62] 


Based on these Go Wo levels, we obtain from the Bethe-Saltpeter 
equation the imaginary and real parts of the dielectric function of 
bulk A-Ti 02 for polarization perpendicular (ordinary) and parallel 
(extraordinary) to the tetragonal axis c shown in Figure [2] These 
are comparable to the dielectric functions obtained from reflection 
spectra polarized perpendicular to the a or c-axis at room temper¬ 
ature by Kramers-Kronig transformations. 66 ^ Note that 86% of the 
experimental reflectivity spectra polarized perpendicular to the a- 
axis is parallel to the c-axis) 66 Furthermore, our dielectric func¬ 
tions agree well with those obtained from BSE calculations within 
the Tamm-Dancoff approximation. 46 In particular, we obtain ex¬ 
cellent agreement both in position and intensity for the first bright 
exciton at ~ 4 eV, which is perpendicular to the c-axis. The lowest 
energy BSE Go Wo transition is at 3.73 eV, about 0.12 eV below the 
PBE Go Wo indirect electronic gap of A-TiC> 2 , as shown in Table [l] 
This is significantly higher than the estimated optical band gap of 
3.42 eV reported in ref. E3 

However, we tend to underestimate the real part of the dielectric 
function, shown in Figure [2jb,d). For example, as reported in Ta¬ 
ble [5] the dielectric constant £00 = = 0) is underestimated by 

about 2 in our BSE calculations. This might be remedied by in¬ 
cluding a greater number of transitions within the BSE calculation. 



hco (eV) fao (eV) 


Figure 2. (a,c) Imaginary and (b,d) real parts of the dielectric function 
of bulk A-TiC >2 for polarization perpendicular (a,b) and parallel (c,d) to 
the A-TiC >2 tetragonal c-axis, Im[£L (&>)], Re[ci(m)], Im[c||(m)], and 
Re[£|| (m)], versus energy (ha> g ), in eV. The BSE spectra from this work 
(red) and from ref. [46] (green) are based on Go Wo eigenvalues. The experi¬ 
mental spectra (blue) are obtained from reflection spectra polarized perpen¬ 
dicular to the (a,b) c-axis or (c,d) c-axis by Kramers-Kronig transformation 
from Ref. [66] 

Table 2. Macroscopic Dielectric Constants Soo Perpendicular ( 1 ) and 
Parallel (||) to the Tetragonal c-axis of A-Ti02 and R-Ti02. 

method xc-functional A-Ti02 R-TiC >2 

1 II 1 II 

7.83 a 9.38 fl 

7.69* 8.91* 

7.61 9.09 

7.55* 9.02* 

5.21 6.09 

5.74* 6.77* 

5.60“ 7.11“ 

5.15 c 6.22 e 

5Jl d 1.33 d 

5.19 e 7.04“ 

5.88* 7.14* 

6.84* 8.43* 

fl Ref. [26] *Ref. [67] c Ref. [46] d Ref. [60] e Ref. [68] f Ref. [66] *Ref. [69] A Ref. ES 


In any case, such computationally demanding calculations are 
beyond the scope of the present work. Overall, the agreement ob¬ 
tained for the BSE dielectric function based on Gq Wq eigenenergies 
demonstrates the robustness of the parameters we will use to calcu¬ 
late the G 0 W 0 PDOS for H 2 0. 

Figure [T] depicts schematically the clean and stoichiometric A- 
TiO2(101) surface. For the clean surface, there are two Ti coordi- 
nately unsaturated sites (Ti CU s) and two bridging O atoms (Ob r ) in 
each unit cell. 

PBE Go Wo places the IP for A-TiO2(101) at 7.15 eV, which is 
0.14 eV below that of R-Ti02(l 10). 26 ^This relative ordering is con¬ 
sistent with, albeit significantly smaller than, that measured with 
XPS for the A-Ti0 2 -Ru0 2 -R-Ti0 2 interface of 0.7 ±0.1 evP 
This ordering also agrees with the 0.47 eV difference in IP cal¬ 
culated using a hybrid quantum-mechanical/molecular mechanical 
embedding technique. 73 In these calculations the IP was obtained 
from the total energy difference upon removal of an electron from 
the neutral A-TiC >2 and R-Ti02 embedded cluster models. Our rela- 



LDA 

7.18 

6.81 

RPA 

PBE 

7.06 

6.60 


HSE 

4.91 

4.83 

BSE-G 0 W 0 

LDA 

4.17 

3.45 

PBE 

4.91 c 

4.76 c 

BSE-DFT 

PBE+A 

5.12 J 

4.98^ 

Experiment 

5.73 / 

5.64 f 
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Figure 3. (a-f) 1ML intact and (g,h) Vi dissociated H 2 O adsorbed on coordinately unsaturated Ti sites (H20@Ti CM5 ). A-Ti02(101)/R-Ti02(110) 21 total 
(grey/black) and H?0 projected (blue/green) DOS computed with (b,c) DFT, (d) scQPGWl, and (e,f,h) Go Wo using the (b,d,e,h) generalized gradient ap¬ 
proximation (PBE) 37 and (c,f) range-separated hybrid (HSEj^for the xc-functional. Filling denotes occupation for A-TiO2(101). Energies are relative to 
the vacuum level, E vac . The measured £cbm 4 -^ (black dashed line), measured and coupled-cluster (CCSD(T)) H 2 O gas phase ionization potentials IP 71 (blue 
dotted line), and for each level of theory the calculated gas phase lbi, 3ai, and lb 2 H 2 O levels 21 (marked in cyan) are provided. Note that the fully symmetric 
2ai H 2 O levels lie below -20 eV. (a,g) Charge transfer of about -OAe accompanying deprotonation is represented by arrows, while intermolecular (gray) and 
interfacial (black) hydrogen bonds are denoted by dotted lines. Ti, O, and H atoms are depicted in silver, red, and white, respectively. 


tive ordering is also consistent with that obtained from KS eigenval¬ 
ues using the B3LYP xc-functional of 8.24^ and 8.6 eV 7 ^for A- 
TiO2(101) and R-TiO2(110), respectively. This qualitative agree¬ 
ment provides further support for the reliability of our Go Wo ap¬ 
proach. 

3.2. H 2 0 Intact and V 2 Dissociated on Coordinately Unsat¬ 
urated Ti Sites. For coverages up to 1ML, H 2 O adsorbs molec- 
ularly on the A-TiO2(101) surface, with O bonding to Tw and 
one H forming an interfacial hydrogen bond with 0/ 7r J 8 11 76 " 78 -*as 
shown in Figure [3ja). On R-TiO2(110), the distance between the 
nearest neighboring Ti CU s sites is shorter, allowing additional inter¬ 
molecular hydrogen bonds to form along the [001] direction] 79 82 
as show in Figure |3ja). 

The QP level alignment relative to the vacuum level E vac for 
1ML of H 2 0 adsorbed intact on A-TiO 2 (101) and R-TiO 2 (110)ED 
are shown in Figure[3pb-f). These are compared to the lb 2 , 3ai, and 
lbi levels’ absolute alignment for gas phase H 2 O . 21 Specifically, 
we analyze the dependence of the F^O PDOS on the methodology: 
PBE DFT, HSE DFT, PBE scQPGWl, PBE G 0 W 0 , and HSE G 0 W 0 . 

As was previously found for the H20-R-Ti02(l 10) interface, the 
IP for H2O-A-TiO2(101) is ordered according to the method’s de¬ 
scription of the screening, £^J 21 As shown in Table [3] the IP is 
ordered as PBE G 0 W 0 (6.3 eV) » PBE scQPGWl (6.3 eV) ~ PBE 
DFT (6.4 eV) < HSE G 0 W 0 (6.9 eV) < HSE DFT (7.2 eV). 

Note that the CBM and VBM relative to E vac from PBE 
scQPGWl and PBE Go Wo are essentially the same for H 2 O-A- 
TiC>2(101), but are significantly lower for H2O-R-TiO2(110). This 
is because the dielectric constant employed in both single-shot PBE 
scQPGWl and PBE Go Wo are those obtained from PBE DFT, 
whereas when the QP self energies are “mixed” with the DFT xc- 
potential in each cycle, as for H20-R-Ti02(l 10), the scQPGWl di¬ 
electric constant is significantly reduced relative to PBE DFT. This 
demonstrates that without mixing of the self energy, for the QP 


Table 3. Ionization Potentials IP in eV of 1ML H 2 0@Ti CMS on A- 
TiO 2 (101) and R-TiO 2 (110). 


method 

xc-functional 

H20@Ti CM5 



A-TiO 2 (101) 

R-TiO 2 (110) 

DFT 

PBE 

6.4 

6.2 a 

HSE 

7.2 

7.3° 

scQPGWl 

PBE 

6.3 

6.6 a 

GqWq 

PBE 

6.3 

6.0 a 

HSE 

6.9 

6.5 a 


a Ref. I2TI 


PDOS the PBE scQPGWl procedure provides no advantage over 
PBE Go Wo, as predicted in ref. EH 

Generally, the highest H 2 O PDOS peaks, £p%$ S , follow the same 
ordering as the IPs. This suggests that is pinned to the VBM 

of the H2O-A-TiO2(101) interface. This is also the case for 1ML 
intact H20@Ti CU s on R-Ti02(l lO).^ 1 However, this ordering of the 
IPs is completely different from that found for gas phase H 2 O. 1 — 'In 
this case, the IP is the energy needed to remove one electron from 
the H 2 O lbi level. Here, the IPs increase with decreasing screening 

121 1 

within the methodology until £00 ~ l. 1 — 1 

However, for 1ML intact H20@Ti CU s, the relative alignment 
of the A-TiO2(101) and R-TiO2(110) VBMs differs qualitatively 
with methodology. The IPs obtained from PBE DFT and PBE 
GqWq are higher (~ 0.2 and 0.3 eV) for A-TiO2(101) than for R- 
TiC>2(110). However, the IP from HSE DFT is lower (—0.1 eV) 
for A-TiC>2(101) than R-TiO2(110), while the opposite is true for 
HSE Go Wo (~ 0.4 eV). Thus, independently of the xc-functional 
employed, Go Wo yields higher IPs for 1ML intact H20@Ti CU s on 
A-TiC>2(101) than on R-TiO2(110). This is contrary to our find¬ 
ings for the clean A-TiO2(101) and R-TiC>2(110) surface, and sug- 
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Figure 4. (a) 1ML intact and (b) Vi dissociated H 2 O adsorbed on coordinately unsaturated Ti sites (H20@Ti CM5 ). Total (grey/black) and H 2 O projected 
(blue/green) Go Wo DOS on anatase TiO2(101)/rutile TiO2(110) surfaces and selected orbitals are shown schematically below/above. Energies are relative to 
the valence band maximum, £vbm- Ti, O, and H atoms are depicted in silver, red, and white, respectively. 


gests that H 2 O adsorption inverts the relative positions of the A- 
TiO 2 (101) and R-TiO 2 (110) VBMs. 

Although we find the position of the lower edge of the va¬ 
lence band is only weakly affected by adsorbing H 2 0 on either 
A-TiO 2 (101) or R-TiO 2 (110), the VBM is shifted up by about 1 
eV in both cases. This is consistent with the experimentally ob¬ 
served change in work function for the liquid water-R-TiO 2 (110) 
interface.EH&EDEllIl 

The reordering of the HSE DFT and Go Wo IPs for 1ML intact 
H 2 0 on A-TiO 2 (101) and R-TiO 2 (110) may be attributed to the 
greater difference between the constant screening built into HSE 
DFT^and the screening of rutile compared to anatase. Essentially, 
the fraction of the Hartree-Fock exact-exchange which is incorpo¬ 
rated within HSE, a = 0.25, acts as an effective inverse dielectric 
constant within the system, a ~ s ~ 0 o 1 . 47 As a result, for materials 
with £00 * 4, HSE DFT and Go Wo should provide similar descrip¬ 
tions of the screening. 4 -^From Table [ 2 ] we see that the RPA, BSE, 
and measured £00 agree qualitatively and are consistently lower and 
closer to the HSE DFT effective dielectric constant of £00 ~ 4 for 
A-Ti0 2 compared to R-Ti0 2 . For this reason, as shown in Fig¬ 
ure [3] the difference between HSE DFT and Go Wo IPs is larger 
for R-Ti0 2 than A-Ti0 2 , resulting in their relative reordering at the 
Go Wo level. This demonstrates the important role played by the 
screening in describing the relative alignment of anatase and rutile 
polymorphs. 

Overall the H 2 0 QP PDOS for 1ML intact H 2 0@Ti C us is similar 
for the A-TiO 2 (101) and R-TiO 2 (110) surfaces. In particular, the 
most strongly bound lb 2 peaks and the upper edges of the H 2 0 
PDOS spectra have similar energies for the two polymorphs over 


all five levels of theory (cf Figure |3jb-f)). 

On A-TiO 2 (101), the 1ML intact H 2 0 QP PDOS generally con¬ 
sists of three distinct peaks, which have clear contributions from 
molecular lb 2 , 3ai and lbi levels ( 1 cf Figure[3jb) and FigureQa)). 
This is in contrast to R-TiO 2 (110), where the H 2 0 QP PDOS con¬ 
sists of many more peaks, with a greater hybridization at 1ML 
compared to 'AML coverage on the R-TiO 2 (110) substrate. 1 — This 
may be attributed to stronger intermolecular interactions on R- 
TiO 2 (110) due to its shorter T\ CU s nearest neighbor separations 
(d[Ti C us -Ticus] » 2.96 A) versus A-TiO 2 (101) (d[Ti CU s - Ti CU s] « 
3.78 A). This leads to intermolecular bonding and antibonding lev- 
els, which may further hybridize with the substrate. 1 — 1 For exam¬ 
ple, as shown in Figure [3] the bottom edge of the 3ai peak for A- 
TiO 2 (101) is higher than that of R-TiO 2 (110). This is because on 
R-TiO 2 (110) the 3ai levels of neighbouring molecules hybridize 
to form intermolecular bonding and antibonding combinations. 1 — 1 
These give rise to separate peaks below and above the bottom edge 
of the R-TiO 2 (110) valence band. As a result, the QP H 2 0 PDOS 
for 1ML intact H 2 0@Ti CM5 on R-Ti0 2 (l 10) has the 3ai intermolec¬ 
ular bonding level below the bottom of the valence band, while for 
A-TiO 2 (101), the 3ai level is completely within the substrate’s va¬ 
lence band. 

Figure [3jg) shows the structures of Vi dissociated H 2 0@Ti CM5 
on A-TiO 2 (101) and R-TiO 2 (110). In both cases, one proton from 
H 2 0@Ti cus is transferred to the adjacent 0^ r . This results in two 
distinct OH groups: HO@Ti CM5 and 0^ r H. This process is accom¬ 
panied by a -0.4e charge transfer from HO@Ti C us to 0/, r H, as de¬ 
picted schematically in Figure [3jg). 

Although the resulting PBE Gq Wq QP DOS shown in Figure|3je) 
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(a) A-TiO 2 -i /4 (101) 
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Figure 5. ’AML H 2 O dissociated on bridging O vacancies (H20@0^ c ) of defective (a/b) A-Ti02-i/ 4 ( 101)/R-TiO2-i/ 4 (l 10)surfaces with V 2 ML of O v h ac . 
Total (grey/black) and H 2 O projected (blue/green) DOS computed with (c,d) DFT and (e,f) Go Wo using the (c,e) generalized gradient approximation (PBE)ED 
and (d,f) range-separated hybrid (HSE) 39 for the xc-functional. Filling denotes occupation for A-TiO2_i/ 4 (101). Horizontal black lines denote Fermi levels s? 
for R-TiO2-i/ 4 (110). Energies are relative to the vacuum level, E vac . Ti, O, and H atoms are depicted in silver, red, and white, respectively. 


and (h) are generally similar, there are important differences which 
are related to the H 2 0@Ti CM s dissociation. For the dissociated H 2 0 
molecule, the lb 2 peak is replaced by separate HO@Ti CMi? and 0^ r H 
peaks at about -12.5 and -13 eV below E vac (see Figure[3jh)), with 
O-H cr character on both A-TiO 2 (101) and R-TiO 2 (110 ) 21 (see 
Figure |4jb)). These peaks are at such similar energies on both A- 
TiO 2 (101) and R-TiO 2 (110) because they are well separated from 
the bottom edge of the Ti0 2 valence band. 

As mentioned above, the three distinct peaks in the H 2 0 PDOS 
on both A-TiO 2 (101) and R-Ti0 2 (l 10) are associated with the lb 2 , 
3ai, and lbi gas phase H 2 0 levels. This is clearly seen by com¬ 
paring the molecular components of the orbitals depicted for 1ML 
intact and l A dissociated H 2 0@Ti CU s on both A-TiO 2 (101) and R- 
TiO 2 (110) in Figure [4] with the gas phase H 2 0 levels depicted 
in Figure [3] There is significantly greater hybridization between 
the molecular levels on R-TiO 2 (110) compared to A-TiO 2 (101). 
Specifically, on R-TiO 2 (110) there are obvious bonding and an¬ 
tibonding combinations of the lb 2 levels and 3ai levels between 
neighbouring H 2 0 molecules. 1 —’Such intermolecular hybridization 
does not occur for A-TiO 2 (101), as the molecules are too far apart. 

More importantly, £p^ S is shifted to higher energy upon disso¬ 
ciation, with a greater shift for A-TiO 2 (101) versus R-TiO 2 (110). 
To explain these differences, and their potential impact on the inter¬ 
faces’ photocatalytic activity, one should compare the level align¬ 
ment relative to the VBM. In so doing, one can directly compare 
the relative propensity for photogenerated hole transfer from the 
substrate’s VBM to the molecular HOMO for A-TiO 2 (101) and R- 
Ti0 2 (l 10). 

In Figure [4] we provide the level alignment relative to the VBM 
for (a) intact and (b) ! /2 dissociated H 2 0 on A-TiO 2 (101) and 
R-TiO 2 (110). The level alignment shown in Figure [4] suggests 
that (1) hole trapping is more favored on A-TiO 2 (101) than R- 
TiO 2 (110) and (2) HO@Ti CM s is more photocatalytically active 
than intact H 2 0@Ti CU s- This is based on the following observa¬ 
tions: (1) £p^ak S is about 0.5 eV higher in energy for A-TiO 2 (101) 


than R-TiO 2 (110); (2) is about 0.1 eV closer to the VBM 

for HO@Ti CM s than for intact H 2 0@Ti C us\ (3) the PDOS for 
HO@Ti CMiS at £pe C ® S is an order of magnitude greater than for intact 
H 2 0@Ti cus\ and (4) the HOMO is more localized on the molecule 
for HO@Ti CM s than for intact H 2 0@Ti C us- 

These conclusions are reinforced by analyzing the HOMOs at 
T shown in Figure [4] Here, one clearly sees that the HOMOs 
have greater weight on the molecule for HO@Ti CM5 than intact 
H 2 0@Ti CM s. This should promote hole trapping on HO@Ti CM £. Al¬ 
though there is only a small (0.1 eV) energy difference between the 
HOMO for ¥i dissociated and intact H 2 0@Ti C us, the latter level is 
not photocatalytically relevant for hole trapping on the molecule. 
This is because it is a lone-pair orbital that datively bonds to Ti CU s- 
For this reason, if an electron were extracted from this level, one 
would instead expect the hole to remain on the surface, and H 2 0 to 
desorb from Ti CU s- This agrees with previous studies of the liquid 
H 2 O-A-TiO 2 (101) interface, which found that localizing the hole 
on intact H 2 0 is inherently unstable, and leads to deprotonation . 17 
Instead, the hole localizes on 3-fold coordinated surface O (O 3 f Q id) 

rm 

atoms “ 

In contrast to the intact H 2 0@Ti CMl? HOMOs, the HOMOs for 
HO@Ti cus on A-TiO 2 (101) and R-TiO 2 (110) are the photocatalyt¬ 
ically active levels for hole-trapping. Indeed, they have the same 
character as the hole trapping levels reported in the literature for 
A-TiO 2 (101 ) 17 and R-TiO 2 (110 ). 22 In particular, they have both 
O 3 fold 2pJ 84 and OH 2p character. While in the case of HO@Ti CM5 , 
this orbital is doubly occupied, in the trapped hole structures of 
refs. [IT] and [22] the OH groups are bent towards the surface, with 
the hole shared between O 3 f 0 id 2p n and OH 2p orbitals. 

This clearly demonstrates that a HOMO initially below the 
VBM can, upon light absorption and subsequent nuclear relaxation, 
evolve into a hole trapping level of the interface. This justifies our 
use of ground state level alignment for comparing photocatalytic 
activity among H 2 0-Ti0 2 interfaces. 

Although hole trapping has been documented for both A- 
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TiO2(101) 17 and R-TiO2(110), ^^ the more favorable ground 
state HOMO level alignment for A-TiO2(101) may explain why the 
anatase polym orph show s higher photocatalytic activity than the ru- 

tile polymorphP^SIl 

3.3. H 2 0 Dissociated on Bridging O Vacancies. For R- 

Ti 02 (l 10), the most stable O vacancies are at surface O br sites, i.e., 
O v b a r c . These sites mediate H 2 O dissociation on R-TiO2(110). 31 88 
For A-TiO2(101), the most stable O vacancies are subsurface. 6 89 
However, after H 2 O adsorption, these subsurface vacancies migrate 
to the surface and are filled by H 2 O, i.e., H 20 @ 0 ^ c , which subse¬ 
quently dissociates to form 2HO^ r J- 89 90 This results in a structure 
equivalent to H on a stoichiometric A-TiO2(101) surface. 89 91 For 
this reason, we consider a ViML coverage of H 2 O adsorbed disso- 
ciatively on O v b a r c sites (H20@0^ c ) of a defective A-TiO2-i/ 4 (101) 
or R-TiO2_i/ 4 (110) 21 surface consisting of V 2 ML of 0^ c , shown 
schematically in Figure[5] This is equivalent to 1ML of H adsorbed 
on O br (H@0^, r ) of a stoichiometric A-TiO2(101) or R-TiO2(110) 
surface. 

These hydroxylated structures have occupied Ti 3 d levels which 
are associated with reduced Ti 3+ atoms. The excess electrons intro¬ 
duce n- type doping. These occupied Ti 3+ 3 d levels give rise to the 
charge density iust below the Fermi level, s?, in the DOS shown in 
Figure |ijc-f). 83 

The difference in spatial and energetic localization of the Ti 3+ 
3 d levels between O defective A-TiO2(101) and R-TiO2(110) has 
been recently probed via STM. 92 ^ For O^ c @A-TiO2(101) at 6 K, 
the excess electrons are strictly localized next to 0^ c ,® while 
for O^ c @R-TiO2(110) at 78 K, the excess electrons are not con¬ 
fined next to O v b a r c ^ 93 Instead, the excess electrons in O v b a r c @R- 
TiO2(110) may occupy 3d levels of surface Ti CU s or subsurface Ti 
atoms. 

We find for YtML U 2 0@0 v b a r c on both A-TiO 2 _i/ 4 (101) and R- 
TiO2_i/ 4 (110), the highest energy occupied Ti 3+ 3d levels 84 are 
mostly on surface Ti atoms, as shown in Figure [5ja,b). These pre¬ 
dominantly Ti 3d x 2_ y 2 levels 84 are bonding along the [010] and 
[001] directions for A-TiO2_i/ 4 (101) and R-TiO2_i/ 4 (110), respec¬ 
tively. Furthermore, for Yi 2 0@O v br c on A-TiO2_i/ 4 (101), the level 
occupies HO^’s nearest neighbor Ti atoms. For H20@0^ c on R- 
TiO2_i/ 4 (110), this level also has weight on the next next nearest 
neighbour Ti CU s atoms. Additionally, there are higher energy occu¬ 
pied Ti 3d levels on subsurface Ti atoms. 

In PBE DFT, the occupied Ti 3d levels form a shoulder 
at the bottom edge of the conduction band for Yi 2 0@O v br c on 
A-TiO 2 -i/ 4 ( 101 ), whereas on R-TiO 2 _i/ 4 ( 110 ) they do not even 
form a shoulder, as shown in Figure [5jc). The degree of 
energetic localization of the Ti 3+ 3d levels, and their energy 
£ Ti 3 + below £p, increases with the level of theory from PBE 
DFT < HSE DFT (s Ti 3+ ~ 0.6,0.4 eV) < PBE G 0 W 0 (s Ti 3+ ~ 
0.7,0.6 eV) < HSE G 0 W 0 (s Ti 3+ ~ 1.0,0.9 eV) for H 2 0@0^ c on 
A-Ti 02 -i/ 4 ( 101 )/R-Ti 02 _i/ 4 ( 110 ), and is generally higher (0.1 eV) 
for A-TiO 2 _i/ 4 ( 101 ) than R-Ti 02 _i/ 4 ( 1 10), as shown in Table[4] This 
is consistent with the s Ti 3 + ~ 1 eV measured for O defective A- 
TiO 2 ( 101 ) and R-TiO 2 ( 110 ) and HO^ r @R-TiO 2 ( 110 ) by scanning 
tunneling spectroscopy (STS)J® 94 photoemission electron spec- 
troscopv (PES), 1 — 1 — 1 and two photon photoemission spectroscopy 
(2PP). 83 96 However, a full treatment of Ti 3d defect levels, e.g., 
due to interstitial Ti atoms, also requires the inclusion of electron- 
phonon interactions. 1 — 1 — 1 

Overall, relative to E vac , the levels of A-TiO2_i/ 4 (101) are consis¬ 
tently about 0.6 eV lower in energy than those of R-TiO2_i/ 4 (110), 
for PBE DFT, HSE DFT, PBE G 0 Wo, and HSE G 0 Wo, as shown 
in Figure |5Jc-f). However, the H20@0^ c lb 2 levels are at similar 
energies (within 0.2 eV) on A-TiO2_i/ 4 (101) and R-TiO2_i/ 4 (110), 


Table 4. Occupied Ti 3+ 3d Level Energies e Ti 3 + 

in eV Below the Fermi 

Level ep for V 2 ML Dissociated H 2 0@0)'" r on 

A-TiO 2 _i/ 4 (101) and R- 

TiO 2 _i/ 4 (110). 



method 

xc-functional 

h 2 o@o^ c 



A-TiO 2 -i/ 4 (101) 

R-TiO 2 _i/ 4 (110) 

DFT 

PBE 

0.2 

0A a 

HSE 

0.6 

0A a 

G 0 W 0 

PBE 

0.7 

0.6 a 

HSE 

1.0 

0.9 a 


STS 

1 . 0 ± 0 . 1 * 

0.7 ±0.1^ 

0.9 C 


PES 

l.l d 

p p 

OO 


2PP 


0.9 s 

a Ref. 1271 b 

Ref. \92\ c Ref. [93] d Rof. [95] 

e Ref. [32] ^Refs. [83landl96l 


for PBE DFT, HSE DFT, PBE Go Wo, and HSE Go Wo, as shown in 
Figure |5Jc-f)- 

Focusing on the IP from PBE Go Wo shown in Table [5] they 
are ordered: R-TiO 2 (110) (7.29 eV) » A-TiO 2 (101) (7.15 eV) 

> H 2 0@Ti a „ on A-TiO 2 (101) (6.25 eV) « H 2 0@0^ c on A- 
TiO 2 _i/ 4 (101) (6.19 eV) » H 2 0@Ti,,, s on R-TiO 2 (110) (6.03 eV) 

> H 2 0@0^ e on A-TiO 2 _i/ 4 (101) (5.37 eV). On the one hand, there 
are no significant differences in IP between bare and H 20 @Ti C us 
covered A-TiO 2 ( 101 ) and R-TiO 2 ( 110 ). On the other hand, 
for H 2 O dissociatively adsorbed on O v br c sites of defective A- 
TiO 2 _i/ 4 ( 101 ) and R-TiO 2 _i/ 4 ( 110 ), the IPs are significantly differ¬ 
ent, despite having quite similar HO^ r concentrations per unit area. 
The origin of this difference might be related to differences in the 

[02I 

structure’s relative stability or surface dipole . 1 — 1 


Table 5. Ionization Potentials IP in eV from PBE GoW q for A-TiO2(101) 
and R-TiO 2 (110) 


coverage 

surface 

IP (eV) 

clean 

A-TiO 2 (101) 

7.15 

R-TiO 2 (110) 

1.29 a 

1ML H 2 0@Ti CM , 

A-TiO 2 (101) 

R-TiO 2 (110) 

6.25 

6.03^ 

1ML H 2 0@0^ c 

A-TiO 2 _i/ 4 (101) 

R-TiO 2 _i/ 4 (110) 

6.19 

5.31 b 

a Ref.|26l^Ref.l2ll 



Similarly, 

for dissociatively adsorbed H 20 @ 0 ^ c on A- 


TiO 2 _i/ 4 ( 101 ) is about 0.4 eV below that on R-TiO 2 _i/ 4 ( 110 ). Since 
£ peak S thus closer to the standard hydrogen electrode (SHE) for 
H 20 @ 0 ^ c on R-TiO 2 _i/ 4 ( 110 ) than A-TiO 2 -i/ 4 ( 101 ), one would 
expect the former structure to require a smaller overpotential and be 

more active than the latter within an electrochemical cell.® How- 

PDOS 

ever, for photocatalysis, the alignment of s peak relative to eyBM 
the more relevant quantity. As we shall see in the next section, the 
relative electrochemical and photocatalytic activities of these two 
structures are reversed. 

3.4. Coverage and Dissociation Dependence of H 2 0 
Spectra for Stoichiometric and Defective Surfaces. To sys¬ 
tematically investigate the similarities and differences between 
A-TiO 2 ( 101 ) and R-TiO 2 ( 110 ) surfaces, we consider a variety 
of coverages of intact and dissociated H 2 0 on stoichiometric 
A-TiO 2 (101) [Figure |6ja)] and defective A-TiO 2 _ v (101) [Fig¬ 
ure ( 6 fb,c)], as done previously for the rutile surface . 21 These con- 
figurations are consistent with previous results for H 9 O on A- 

Ti q 2 j8inr76it2gE) 

The adsorption energies shown in Tableland Figure |7Ja,b) for 
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a) A-Ti0 2 (101) - Ti cus (b) A-Ti0 2 _, /8 (101) - VML O If 



Figure 6. Schematics of H 2 O adsorbed intact (I) or dissociated (D) on 
(a) coordinately unsaturated Ti sites (Ti cus ) of stoichiometric A-TiO2(101) 
and (b) !4ML or (c) VML of bridging O vacancies (O v b a r c ) on defective A- 
TiO2- x (101), where x = Vs and V 4 , respectively. Coverage is half the number 
of H 2 O formula units per (101)1x1 unit area of the clean stoichiometric or 
defective surface. Dissociation is the fraction of H 2 O molecules which are 
dissociated. Charge transfer of about -OAe accompanying deprotonation is 
represented by arrows, while intermolecular (gray) and interfacial (black) 
hydrogen bonds are denoted by dotted lines. 


H 2 O on A-Ti 02 , A-Ti 02 _i/ 4 , and A-Ti 02 _i/ 8 (101) are generally 
similar to those on R-T 1 O 2 , R-Ti02-i/ 4 , and R-Ti02-y 8 (110)“ 
respectively. On both A-TiO2(101) and R-TiO2(110), intact H 2 O 
adsorption is more stable than dissociative adsorption from Vi and 
W2 ML coverages. The adsorption energies for H20@Ti C us on A- 
TiO 2 ( 101 ) follow the same trend as on R-Ti 02 (l 10), but are some¬ 
what stronger on A-TiO2(101), with the greatest differences seen 
for dissociatively adsorbed H 2 O. Since the photocatalytically ac¬ 
tive species HO br @Ti CU s is more stable on A-TiO2(101) than R- 
TiO2(110), this also suggests that A-TiO2(101) should be more 
photocatalytically active than R-TiO2(110). This is because one 


Table 6. Adsorption Energies E a ^ s , Highest PDOS Peaks and 

Average PDOS HOMO Energies ^qaio * n H 2 O on Ticus of Sto¬ 
ichiometric A-TiO2(101) and R-TiO2(110) and 0™ c of Defective A- 
TiO 2 -jc(101) and R-TiO 2 _,(110) with x = Vs or Va. 


coverage A-TiO2-jc(101) R-TiC^-j^llO) 


ML 

X 

Eads 

PDOS 

s HOMO 

PDOS 

£ peak 

Eads 

PDOS 

s HOMO 

PDOS 

£ peak 

V2 

I 

0 

-0.38 

-1.42 

-0.7 

-0.33 

-1.44 

-1.2 

D 

0 

-0.18 

-0.49 

-0.7 

-0.13 

-1.03 

-0.8 


I 

0 

-0.37 

-1.07 

-0.6 

-0.41 

-1.28 

- 1 . 1 * 

1 

V2D 

0 

-0.29 

-0.43 

-0.5 

-0.23 

-0.79 

-1.0* 


D 

0 

-0.26 

-0.45 
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expects there to be more HO^ r @Ti CMi? on A-TiO2(101) than R- 
TiO 2 (110). 

On the defective surfaces, intact H 2 O adsorption is also more sta¬ 
ble than dissociative adsorption on A-TiO2_i/ 4 (101)/R-TiO2_i/ 4 (l 10) 
and A-TiO 2 -i/ 8 (101)/R-TiO 2 -i/ 8 (l 10). Dissociative H 2 0@0^ c ad¬ 
sorption is generally stronger on A-Ti02-y 4 (101)/A-Ti02-y 8 (101) 
than R-Ti02_i/ 4 (110)/R-Ti02_i/ 8 (110) surfaces, except for WML 

1 / 3 DH 2 0@0^ c . 

However, the adsorption energies shown in Figure 0b) 
are strongly dependent on the stability of the defective A- 
TiO 2 -i/ 4 (101)/R-TiO 2 -i/ 4 (l 10) and A-Ti02-i/ 8 (101)/R-Ti0 2 -i/ 8 (110) 
structures with surface O v b a r c . Since surface O v br c are less stable 
than subsurface O vacancies for A-TiO2(101), the adsorption en¬ 
ergies on A-TiO2-*(101) provided in Figure [Tpb) are somewhat 
overestimated. 

Figure[7Jc,d) shows the PBE Go Wo H 2 O PDOS relative to £vbm 
as a function of coverage and dissociation for the structures shown 
in Figure [6] Overall the PDOS on A-TiO2(101) and R-TiO2(110) 
are in surprisingly close agreement, both in shape and energy. For 
VML of H 2 O, peaks related to the H 2 O lb 2 , 3ai and lbi levels (cf 
Figure [3jb» are clearly evident. When the coverage is increased 
to more than 1ML, there are larger differences between the H 2 O 
PDOS on A-TiO2(101) and R-TiO2(110). This may be attributed 
to the different intermolecular and interfacial interactions induced 
by the different hydrogen bonding networks. For IV2 ML H 2 O on 
A-TiO2(101), the peak associated with the H 2 O lb 2 level, which is 
located at -8 eV, is more delocalized than rutile. This is because 
there are more interfacial interactions between H 2 O and A-Ti 02 
(101). For lVML H 2 O on 0^ c ( * 1 /sD), the water lb 2 level splits into 
two peaks, as H 2 O and HO^ r form two lines of hydrogen bonding 
networks. We also find the bottom of the valance band for H 2 O 
on A-Ti02, A-Ti02_i/ 4 , and A-Ti02_i/ 8 (101) is higher than that of 
R-Ti02, R-Ti02-i/ 4 , and R-Ti02_y 8 (HO). This is attributable to the 
higher YBM of the clean A-TiC>2(101) surface, as depicted by the 
gray regions in Figure [71 

* I - 1 PDO S 

For all spectra shown in Figure r7fc, d), %omo * s higher on A- 
Ti 02 , A-Ti 02 _i/ 4 , and A-Ti 02 _i/ 8 (101) than R-Ti 02 , R-Ti 02 _i/ 4 , 
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Figure 7. Structure and coverage dependence of (a,b) adsorption energy E a d s and (Cjd) Go Wo PDOS for H?0 adsorbed intact (I) or dissociated (D) on (a,c) 
coordinately unsaturated Ti sites (Ti CM5 ) of stoichiometric A-TiO2(101) (blue, Figurepla)) and R-Ti02(l 10)(green) and (b,d) bridging O vacancies ( O v b a r c ) 
of defective A-TiO2- x (101) (blue) and R-TiO2- x (110) 21 (green), with x = Vs (thin lines, Figure [6jb)) or 14 (thick lines, Figure [SJc)). (a,b) RPBE E ac [ s on 
A-TiO2-;c(101) (O) and R-TiO2-;c(110)^(O) surfaces (x = 0, Vs, 14) for (white) low (14 and 14ML), (turquoise) medium ( 3 A and 1ML), and (blue) high (114 
and 114ML) coverage. (c,d) Energies are relative to the VBM (^vbm)- The clean surface DOS of (c) A-Ti02(101)/R-Ti02(110) (dark/light gray regions) are 
shown for comparison. 
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Figure 8. Average energy of Go Wo PDOS HOMO in eV of 

H20@Ti cus on stoichiometric A-TiO2(101) versus R-TiO2(110) and of 
H20@0 ^ c on defective A-TiO2-jt(101) versus R-TiO2- x (110) for x = Vs 
or 14. H 2 O total coverage in ML and fraction intact (I) or dissociated (D) 
are provided. A linear fit (red dashed line) with a standard deviation of ±0.1 
eV (gray regions) is compared to the identity line (black solid line). 


and R-Ti 02 -i / 8 (110), respectively, as shown in Table [ 6 ] and Fig- 

I I PDOS ^ 

ure[ 8 | Further, these differences in £ H qmo are l ar g er f° r dissociated 

H 2 O. Since it is the HOMO of HO@Ti C us which can trap a photo¬ 
generated hole, as discussed in section [3^2] the larger differences in 
^homo shown in Figuresj^c) andj^for dissociated H 2 O suggest A- 
TiO 2 ( 101 ) should generally be more photocatalytically active than 
R-Ti 02 (l 10) from low coverage (!4ML H 2 O) to multi-layered H 2 O 
(F/ 2 MLH 2 O). 

For V 2 ML of dissociatively adsorbed H 20 @ 0 ^ c , s P p ®ak relative 
to the VBM for R-Ti0 2 _i/ 4 (110) is below that for A-TiO 2 _i/ 4 (101). 
This suggests HO^ r @0^ c should be more photocatalytically ac¬ 


tive on A-TiO 2 _i/ 4 ( 101 ) compared to R-TiO 2 _i/ 4 ( 110 ). However, as 
shown in the previous section, the reverse is true for their relative 
electrochemical activity, i.e., HO^ r @0^ c on R-TiO 2 _i/ 4 ( 110 ) is ex¬ 
pected to be more electrochemically active than A-TiO 2 _i/ 4 ( 101 ). 
This demonstrates the importance of considering both the absolute 
level alignment relative to E vac , and the level alignment relative to 
£VBM- 

4. CONCLUSIONS 

In heterogeneous catalysis, photocatalytic activity is controlled by 
the level alignment of the adsorbate and substrate levels. For this 
reason it is essential to obtain a quantitative description of the in¬ 
terfacial level alignment to determine and predict catalytic activity. 
This can only be obtained from many-body QP GW calculations, 
which are necessary to correctly describe the anisotropic screening 
of electron-electron interactions at the catalyst’s interface. 

Previously, we have shown that HSE GqWq reliably describes 
the interfacial level alignment relative to the VBM for highfy hy¬ 
bridized and localized molecular levels of H 2 O 21 and CH 3 OH 26 on 
R-TiO 2 ( 110 ). Here, we have shown that HSE Go Wo also provides 
a quantitative description of the occupied Ti 3+ 3 d level’s alignment 
relative to the Fermi level on both reduced anatase and rutile poly¬ 
morphs. These are the levels from which electrons are typically 
excited in 2PP experiments .® 5 83 96 Since HSE DFT fails in both 
cases, these results clearly demonstrate the important role played 
by anisotropic screening of the electron-electron interaction in de¬ 
scribing the alignment of these molecular and defect levels. 

In this study we have performed an in-depth comparison of 
the QP GW level alignment for H 2 O-A-TiO 2 ( 101 ) and H 2 O-R- 
TiO 2 ( 110 ) interfaces for a range of chemically significant struc¬ 
tures. We have considered the limits of low and high H 2 O coverage, 
intact to fully dissociated H 2 O, and stoichiometric to O defective 
surface composition. Using the HOMO-VBM level alignment for 
these systems prior to irradiation £^ofvK)’ we ^ ave esta blished the 
following trends in their relative photocatalytic activity for H 2 O 
photooxidation. (1) There is a strong linear correlation between 
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^homo on A-TiO2-jc(101) and R-TiC^-jcCllO). (2) We consistently 
find H 2 0 ’s ^homo c l° ser t0 £ vbm for A-Ti02 than R-Ti02* (3) 

PDOS 

These differences in £homo are g reater for dissociated H 2 O, and in- 
crease as £hqmo approaches sy bm- (4) Overall, %omo approaches 
£vbm with H 2 O dissociation. Altogether, this suggests HO@Ti CM5 
is more photocatalytically active than intact H 20 @Ti CM5 and hole 
trapping is more favorable on A-TiO2(101) than R-Ti02(l 10). This 
may explain why the anatase polymorph is generally more photo¬ 
catalytically active than rutile for H 2 O photooxidation. 

We have clearly demonstrated that the ground state interfacial 
level alignment is a key factor in understanding the photocatalytic 
activity of Ti02. Moreover, in general, knowledge of an interface’s 
ground state electronic structure can be used to establish trends for 
predicting photocatalytic activity. 
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